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ABSTRACT 

We present HST/ GHRS echelle observations of multiple interstellar lines of C I, 
Mg I, Cr II, and Zn II towards both stars in the [i Cru binary system. Despite large 
differences in the profiles of the neutral species, no significant variations between the 
stars are seen in the Cr II and Zn II line profiles. In particular, the Zn II absorption 
observed at -8.6 km/sec towards \i Cru is constant despite greatly enhanced columns 
of the neutral species at this velocity towards /x 1 Cru. An analysis of the fine-structure 
excitation of C I in this cloud implies that the density is njj < 250 cm~ 3 . From the 
lack of variation in the (optical) Na I D2 line profiles towards n 1 and fi 2 Cru in spectra 
taken 21 months apart, we can place a lower limit to the size of the structures of ~10 
AU. These results are discussed in the context of recent radio and optical studies of 
apparently pervasive high density small-scale interstellar structure. 



Subject headings: ISM: structure - ISM: clouds - stars: individual (HD 112092, HD 
112091) 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by the Association for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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1. Introduction 

There is growing evidence from recent radio and optical observations that the diffuse 
interstellar medium (ISM) exhibits significant subparsec-scale variations down to limits of a 
few AU (Frail et al. 1994; Davis et al. 1996; Meyer & Blades 1996; Watson & Meyer 1996). 
Utilizing multi-epoch observations of 21 cm absorption towards high- velocity pulsars, Frail et al. 
find pervasive variations in the H I opacity on scales of 5 to 100 AU with implied densities of 
nn ~ 10 4 — 10 5 cm -3 . Probing larger scales (500-30,000 AU) through optical observations of the 
interstellar Na I absorption towards a number of resolved binary star systems, Watson & Meyer 
(1996) also found a ubiquitous pattern of small interstellar structures and inferred densities of 
n H ^ 10 3 cm~ 3 in some structures. Due to their large overpressures with respect to the intercloud 
medium, such structures cannot be accommodated in any abundance in the standard McKee & 
Ostriker (1977) model of the ISM. Possible explanations for these structures range from small-scale 
filamentary geometries (Heiles 1997), to fractal geometries driven by turbulence (Elmegreen 1997), 
to a separate population of small, dense self-gravitating clouds (Walker & Wardle 1998). 

An important point discussed by Heiles (1997) is that the high volume densities of the 
small-scale structures have not been measured directly - they have been inferred by dividing the 
observed column density differences by the sampled transverse length scales. Furthermore, in the 
case of the Na I observations, since Na I is not the dominant ion in H I clouds it may not be a 
reliable tracer of H I column density differences, especially at low column densities (Welty, Hobbs, 
& Kulkarni 1994). In this Letter, we present HST GHRS echelle observations of the interstellar 
absorption for C I, Mg I, Cr II, and Zn II towards both members of the binary \x Cru (Meyer 
& Blades 1996) in the first effort to apply UV interstellar diagnostics towards an understanding 
of the physical conditions in the structures giving rise to the observed small-scale variations. In 
addition, we present new high-resolution Na I observations of this binary that afford a temporal 
probe of the ISM structure at a scale similar to that sampled by Frail et al. (1994). 

2. Observations 

2.1. HST/GUKS Observations 

The [i Cru binary system consists of a B2IV-V and a B5Ve star at a distance of 170 pc 
with a separation of 38". 8 (corresponding to a projected linear separation of 6600 AU or 0.03 
pc). Both stars were observed using the HST Goddard High Resolution Spectrograph (GHRS). 
Four wavelength intervals were covered using the high resolution echelle gratings (ECH-A near 
127lA, and ECH-B near 2026, 2062, and 2852A) to observe lines of C I, Mg I, Cr II, and Zn II 
at a resolution of ~3.5 km/sec. The data were taken in FP-SPLIT mode, the calibrated data 
files as delivered by STScI were coadded in the usual fashion using the STSDAS tasks poffsets, 
dopoff , and specalign. For the rapidly rotating star fi 2 Cru the continuum was fitted using 
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low-order polynomials. For /i 1 Cru continuum fitting was more difficult due to the low v ■ sin(i) 
(48 km/sec). Here smaller regions of the spectrum were fitted using higher-order polynomials; the 
resulting profile for each line was then compared to the other line of the same species towards 
/i 1 Cru (and the same line towards [i 2 Cru) to be sure that the continuum was properly placed 
and that the entire interstellar line was excluded from the fit. The resulting profiles are shown in 
Figure 1; note that the variation is largely confined to the lines of the neutral species, with no 
significant variations in the Cr II or Zn II profiles. 

2.2. Anglo- Australian Telescope/UHRF Observations 

The interstellar Na I D2 absorption was re-observed with the Ultra-High Resolution Facility 
(UHRF) on the Anglo-Australian Telescope (AAT) towards both members of the fi Cru binary 
star system in March 1997, 21 months after the first observations reported by Meyer & Blades 
(1996). A Tek 1024x1024 CCD was used for these observations; the data was then processed 
using a combination of Figaro^ and IRAFf] routines. The spectrograph configuration for these 
observations is the same as that described in Meyer & Blades (1996); the resulting resolution was 
~ 10 6 . 

According to the Hipparcos catalog, the proper motion of these stars is ~ 0.034"/year (ESA 
1997), corresponding to a projected change in position of ~ 10 AU over the 21 month period 
between observations. As shown in Figure 2, the line profiles are virtually identical; in addition, 
the column densities, b-values, and relative velocities of the various components derived by profile 
fittingF] are well within the calculated errors. The striking lack of change in the Na I D2 profile 
implies a minimum cloudlet size of a few tens of AU. This result is not necessarily inconsistent 
with the ~ 5 — 100 AU structures found by Frail et al. (1994), but suggests that further monitoring 
of this sight-line would be worthwhile. 

3. Discussion 

The narrow, hyperfine-split absorption observed in Na I at -8.6 km/sec towards /x 1 Cru 
provides an opportunity to derive the physical conditions in a relatively isolated component. 
In particular, the components at 4.4 and 5.0 km/sec (which appear to be strongly enhanced 
towards [i 2 Cru) are severely blended at the resolution of the GHRS echelle gratings, and thus 

2 Figaro is distributed as part of the Starlink Project, which is run by CCLRC on behalf of PPARC. 

3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 

4 The fits6p program was used to fit the data, this program is a somewhat modified descendant of the program 
described by Welty, Hobbs, & York (1991) 
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any interpretation of these components is more ambiguous. Profile fitting was used to derive the 
column densities listed in Table 1 for the -8.6 km/sec component. The most remarkable aspect of 
these numbers is the enhancement over typical ISM values of the columns of the neutral species 
relative to the observed Zn II towards Cru, despite the similarity of the Zn II columns. Since 
Zn II is the dominant ion in typical H I clouds, and Zn is relatively undepleted (Sembach et al. 
1995; Roth & Blades 1995), it is a direct measure of the H I column - the lack of variation in 
the Zn II profile implies that the variation in Na I and the other neutrals is occurring without a 
change in the H I column. These observations suggest that we may be seeing a small "core" in the 
cloud at -8.6 km/sec where the neutrals vary, but where any variation in the Zn II (and hence H I) 
column is masked by a broader component with no associated neutral gas. We will now explore 
the physical conditions in this cloud, as well as several potential explanations for the variations, 
and show that small-scale density and temperature fluctuations in a diffuse interstellar cloud are 
the most likely explanation for the observed neutral enhancements towards fi Cru. 

Since fi 2 Cru is a known Be star, one must consider the possibility that the variations in 
absorption line profiles arise in circumstellar gas around one or both members of the system, 
but there is strong evidence that this gas is indeed interstellar in nature. The narrow, hyperfine 
split Na I components observed to vary between the members of the [i Cru system arise in cold 
gas, and it is difficult to imagine an environment near these B stars in which such cold clumps 
would arise. The lack of time-variability in the D2 absorption profiles also argues against a 
circumstellar explanation since an expanding circumstellar shell gas would be likely to shift or 
vary in strength over this period. Finally, the lack of strong C I* (and C I**) absorption in the 
-8.6 km/sec component implies that the radiation field cannot be much larger than ~ 20 times 
the diffuse ultraviolet background (Jenkins & Shaya 1979); calculations show that the observed 
N(C I*)/N(C Itotai) ratio of < 0.25 requires that this structure lies at least 1 pc from the B2V star 
H 1 Cru. 

We can estimate the hydrogen column density, and hence the density for this component using 
the measured uncertainties in the Zn column (towards both members of the system) as a rough 
upper limit to the column of Zn II in the structure giving rise to the Na I variation. We place a 
2cr upper limit on any enhancement at -8.6 km/sec of Zn II towards fi 1 Cru of 4 x 10 10 cm -3 . For 
typical interstellar depletions this corresponds to N(H l)< 2.5 x 10 18 cm~ 3 (Sembach et al. 1995; 
Roth & Blades 1995), which is less than 40% of the value one would estimate from the observed 
Na I column. Given the binary separation of 6600 AU, this implies that the observed variation 
may arise in a structure with a density contrast of nn < 25 cm~ 3 , which is appreciably less than 
that estimated in previous studies of small-scale ISM structure. However, note that 6600 AU 
is only an upper limit to the size of the structure, and a smaller structure would have a higher 
density. Therefore, we now refine this simple argument using the available data for the neutral 
species to directly estimate the density without assuming a characteristic length scale. 

First, we can place an upper limit on the hydrogen density using the fine-structure equilibrium 
of C I. Given the measured ratio of N(C I*)/N(C I toto ;)<0.25, we can place a strict upper limit 
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to the density for cloud temperatures > 20 K of nn < 250 cm~ 3 (Jenkins & Shaya 1979). For a 
more typical interstellar cloud temperature of ~ 100 K the density would be nn < 50 cm -3 , which 
is (again) substantially lower than the densities estimated for small-scale structures in previous 
studies. This limit is independent of any estimate of the cloud size. 

The measured columns of the neutral species can be combined with ionization and 
recombination coefficients to derive estimates of the density in this structure. If we take Zn II as a 
surrogate for the (unobserved) first ions of carbon, sodium, and magnesium, we can estimate the 
required columns by assuming Solar abundance ratios and a variety of depletions (ranging from no 
depletion to typical cold cloud values). Pequignot & Aldrovandi (1986) have calculated ionization 
and recombination coefficients as a function of temperature and electron density for four choices 
of radiation field, including that of Witt & Johnson (1973). For a temperature of 100 K and an 
assumed n e /n# ratio of 2 x 10~ 4 (from the ionization of carbon), we can estimate nn using our 
observed columns (see Table 1). First, to minimize the neutral fraction for C, Na, and Mg we will 
use the total observed Zn II column at -8.6 km/sec in these estimates. Then, for no depletion, we 
derive uh < 205, 120, and 190 cm -3 from C I, Na I, and Mg I respectively. For a typical cold cloud 
depletion pattern (Savage & Sembach 1996), we estimate densities of njj < 260, 240, and 620 cm~ 3 
respectively. Note that the densities estimated from C I and Na I are consistent in both cases, 
but that the density is apparently over-estimated for Mg I if we assume cold cloud depletions. If 
the upper limit of N(Zn II) < 4 x 10 10 cm -3 estimated above is used, densities of nH < 730, 425, 
and 540 cm -3 are found assuming no dust depletion. These density estimates imply cloud sizes 
of ~ 350 - 4,000 AU, which are comparable to the cloud size limits of 10 and 6600 AU set by the 
optical observations. The agreement between the estimated cloud sizes and the observed limits is 
suggestive, but note that there is an implicit assumption that there is little or no ionized hydrogen 
in this structure. If some fraction of the hydrogen is ionized the density in the structure will 
be lower; however, as the ionized fraction increases, the hydrogen column in the structure must 
decrease so as not to exceed the limit on the cloud size set by the binary separation. 

Now in addition to being density sensitive, the ratios of the columns of neutrals to first ions 
is also sensitive to the temperature. For C I, Na I, and Mg I the column density goes roughly 
as ~ x -0 ' 65 (Pequignot & Aldrovandi 1986), thus a difference in temperature of a factor of 100 
would lead to a difference in the columns of the neutrals of a factor of ~ 20, a value consistent 
with the observed upper limit for Na I at this velocity towards fj, 2 Cru. However, the profile fits 
to the Zn II lines towards both stars yield similar b-values for both lines-of-sight suggesting that 
the temperature of the bulk of the material in these clouds is similar and of order a few hundred 
degrees K (given b ~ 0.4-0.5 km/sec), although the presence of a somewhat narrower colder 
"core" towards fi 1 Cru cannot be ruled out. Since the 21 cm opacity scales as T _1 , it seems at 
first glance that temperature fluctuations alone could explain the observed variation in both the 
neutral and the 21 cm lines. But as noted by Heiles (1997), such extreme temperature fluctuations 
are difficult to understand for a number of reasons. This suggests that temperature fluctuations 
are likely to play at least some role in producing the observed neutral variations, although it is 
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unlikely to be the sole mechanism involved. Thus, a plausible model for the absorbing clumps is 
that we are detecting a small, somewhat denser and colder core towards fi 1 Cru that is embedded 
in a larger cloud seen towards both stars. 

There are a number of similarities and differences between the properties of the -8.6 km/sec 
cloudlet towards fi 1 Cru and the structures observed in radio studies. The most direct comparisons 
can be made with the pulsar study of Frail et al. (1994) as it is also a multi-epoch study of 
absorption against a background point-source. The scales sampled in these studies are similar, 
thus the lack of variation observed towards both members of the [i Cru system is somewhat 
surprising given the ubiquity of the variations observed in the pulsar study, but is not inconsistent 
with the inferred cloud sizes. The lack of variation in the observed Zn II profiles between the stars 
suggests that there is little or no variation in the total hydrogen column density towards \i Cru[^, 
which appears to be at odds with the large H I variations observed in the radio studies. The upper 
limit for the density in this structure is significantly different from the values estimated from the 
21 cm observations, which typically imply much larger densities in (perhaps) smaller clouds. One 
explanation is that perhaps this gas is associated with one of the Heiles (1997) filaments or sheets 
observed roughly "face-on", our estimate for the size of the structure would then be an estimate 
for the thickness of the filament /sheet (with the width of the structure being constrained by the 
6600 AU binary separation). A large variation in the Zn II column density could be measured if 
the filament /sheet was along the line-of-sight instead of perpendicular to it. This sort of density 
fluctuation is also consistent with a fractal ISM model (such as those of Elmegreen 1997), since 
density as well as column density inhomogeneities can be a natural outcome of turbulent processes. 

What, if any, fraction of the cloudlets detected as Na I variations are associated with the 
structures responsible for the 21 cm variations towards radio sources is still uncertain. Further 
observations of multiple star systems with variable Na I (in both lightly and more highly reddened 
sight-lines) using HST/STIS and FUSE are required for a fuller understanding of the physical 
conditions in the Na I variable components, and to (potentially) identify components which show 
large H I column density differences similar to the 21 cm studies. If no such components are found, 
then this suggests that the Na I and 21 cm studies are identifying different structures. 

It is a pleasure to acknowledge several valuable conversations with the late Lyman Spitzer, 
with Mark Wardle about extreme scattering events, and with Dan Welty about the ionization 
equilibrium of neutral species. In addition, we would like to thank the anonymous referee for their 
suggestions. 



5 Note that the upper limit to the density from the fine-structure equilibrium of C I and the relatively low columns 
of the observed neutral species imply that Zn I cannot be the dominant ionization stage of Zn in this velocity interval. 
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Fig. 1. — Optical (AAT) and ultraviolet (GHRS) lines observed towards both members of the 
binary system fi Cru, where the solid profile is that observed towards fi 1 Cru, and the dashed is 
that observed towards \j? Cru. Note the relatively large differences between the lines-of-sight in 
the neutral species, compared to the slight differences seen for the first ions. 



Fig. 2. — Comparison of the Na I D2 line profiles observed 21 months apart towards both 
components of the binary system \i Cru. Given the proper motion of fi Cru, the lack of variation 
in these profiles implies that the responsible structures are appreciably larger than 10 AU in size. 



Table 1: Column Densities for the -8.6 km/sec Component towards ji 1 Cru 



Ion 


Vacuum 


Oscillator 


Column 




Wavelength 


Strength 3 - 


Density 




(A) a 




(cm" 2 ) 13 


CI 


1277.2454 


0.09665 


8.0±1.1 x 10 11 


CI* 


1277.2823 


0.07249 


< 2.7 x 10 11 


CI** 


1277.5496 


0.08117 


< 2.7 x 10 11 


Nal 


5891.5833 


0.6550 


6.78±0.08 x 10 1 




5897.5581 


0.3270 


6.94±0.15 x 10 1 


Mgl 


2026.4768 


0.1154 


4.1±0.5 x 10 11 




2852.9642 


1.830 


4.2±1.1 x 10 11 


Cr II 


2062.234 


0.777 c 


< 1.4 x 10 11 


Zn II 


2026.136 


0.5150 


1.4±0.3 x 10 11 




2062.664 


0.2529 


1.3±0.4 x 10 11 



"From Morton (1991) unless noted otherwise. 

^Column densities derived from profile fitting the interstellar absorption corresponding to the -8.6 km/sec component 
towards fj 1 Cru. The errors quoted are la, with 2a upper limits given for the lines of C I*, C I**, and Cr II. 
c Oscillator strength from Bergeson & Lawler (1993). 
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